Renormalization Group Study of the Electron-phonon Interaction in High T c Cuprates 
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We generalize the numerical renormalization group scheme of Ref. [1-3] to study the phonon- 
mediated retarded interactions in the high T c cuprates. We find that three sets of phonon-mediated 
retarded quasiparticle scatterings grow under RG flow. These scatterings share the following com- 
mon features: 1) the initial and final quasiparticle momenta are in the antinodal regions, and 2) the 
scattering amplitudes have a x 2 — y 2 symmetry. All three sets of retarded interaction are driven to 
strong coupling by the magnetic fluctuations around (7r, 7r). After growing strong, these retarded 
interaction will trigger density wave orders with d-wave symmetry. However, due to the d-wave 
form factor they will leave the nodal quasiparticle unaffected. We conclude that the main effect of 
electron-phonon coupling in the cuprates is to promote these density wave orders. 



One of the most prominent puzzles to emerge from 
experimental studies of the high-temperature supercon- 
ducting cuprate materials is the strong momentum-space 
dependence of the properties of the electronic excitations. 
This can be seen most clearly in angle-resolved photoe- 
mission spectroscopy (ARPES) of LSCO: at optimum 
doping, sharp quasiparticle peaks can be distinguished 
near the Fermi surface in both nodal and antinodal di- 
rections. However, when the system is underdoped, the 
antinodal quasiparticle peaks disappear, while the nodal 
quasiparticle peaks remain[4]. Another beautiful demon- 
stration of the robustness of the nodal quasiparticle is 
obtained in the NaCCOC materials. Recent STM stud- 
ies have shown the presence of a commensurate, 4 lat- 
tice constant, checkerboard order which is independent 
of doping [5] (indicating the importance of lattice pin- 
ning). Despite this, ARPES studies of the same system 
reveal nodal quasiparticle peaks [6]. Finally, ARPES ex- 
periments probing the isotope effect in optimally doped 
Bi2212 reveal a substantial isotope shift in the ARPES 
spectra near the antinode, but little shift in the nodal 
spectrum[7]. The broad picture emerging from these ex- 
periments is of nodal quasiparticles which are insensitive 
to doping, disorder, charge and spin order, and lattice 
vibration, while antinodal excitations are sensitive to all 
these perturbations. One of us (DHL) has dubbed this 
the "nodal-antinodal dichotomy" [4] . 

In the cuprates it is widely accepted that in the over- 
doped regime, quasiparticles are well-defined excitations 
all along the normal state Fermi surface. We view the 
effect of decreasing doping as changing the Fermi sur- 
face geometry and increasing the strength of the residual 
quasiparticle interaction. We model this residual interac- 
tion with a momentum independent quasiparticle scatter- 
ing. It is important not to confuse this effective residual 
quasiparticle scattering with the bare local electron cor- 
relation. Previously Honcrkamp et al have performed a 



one-loop renormalization group (RG) study of this resid- 
ual interaction[l]. Here we generalize their method to 
study the phonon-mediated retarded interaction. The 
RG treatment of the retarded interaction is similar to 
that used for one dimensional systems in Ref. [8-10]. Our 
motivation for studying the electron-phonon interaction 
in the cuprates are twofold: 1) we hope to better un- 
derstand the electron-phonon interaction in doped Mott 
insulators, and 2) we hope to gain some insight about the 
origin of the nodal-antinodal dichotomy. 

In a recent paper Devereaux et al pointed out the im- 
portance of the momentum dependence in the electron- 
phonon coupling constant when interpreting ARPES 
data[ll]. In particular they showed that as the ini- 
tial electron momentum is varied along the Fermi sur- 
face, the matrix element that couples the electron to 
the Bi g phonon exhibits a x 2 — y 2 symmetry. Mean- 
while studies including weak to intermediate Hubbard 
interaction (U < 6t) have come to the conclusion that 
the s-symmetry electron-phonon coupling is suppressed 
by electron-electron repulsion [12-1 5], especially for large 
momentum transfer processes [13-15]. 

In this study we follow the RG flow of both the instan- 
taneous and retarded electron-electron scattering ampli- 
tude in the full first Brillouin zone. We compare the 
flow of retarded interaction with different symmetry. We 
start with a tight binding dispersion given by e(k) = 
- 2t [cos ( k x ) +cos ( k y ) ] +4i ' cos ( k x ) cos ( k y ) +4i " [cos 2 ( k x ) + 
cos 2 {k y ) - 1], with H = 0.3i, t" = -0.lt, and fi = -0.7t. 
This set of parameter choice produces a rather realis- 
tic Fermi surface. Following Ref.[l] we choose the bare 
residual quasiparticle scattering amplitude U — St. As in 
Ref. [1] we perform a one-loop Wilsonian RG for the one- 
particle-irreduciblc four-point vertex functions. At each 
energy scale the renormalized four-point function serves 
as an effective interaction for particles with excitation 
energies |e(k)| below the cutoff scale A. The results we 
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report are obtained for a temperature fcgT = 0.04i, and 
the RG flow is integrated between initial cutoff A = At 
and final cutoff 0.2t. With spin-rotational invariance, the 
cutoff-dependent effective interaction is given by 

r 3 

S int = J2 I Yl^k^^VA^MM^H^LUi) 

a.cr' i=l 

l\>\, (k 3 , W 3 ) ?/V (ki , Ul (k 2 , UJ 2 ) . ( 1 ) 

In the above kj) annihilates an electron with quan- 

tum numbers w^k^cr; furthermore o) 4 = W! + w 2 - w 3 
and k 4 = ki + k 2 — k 3 . The contribution to the RG 
flow OaVa is summarized by the Feynman diagrams in 
Fig. (1), where the label "I" or "R" denotes instanta- 
neous or retarded interaction respectively. For the mo- 
ment we discuss the instantaneous interaction, hence all 
interactions are the "I" type. In each diagram, there are 
two internal lines. One represents the Greens function 



G A (k,w) = 



XA(k) 



iuj - e(k) - XA(k)S(k, iu>) ' 
while the other represents 



S A (k,w) 



XA (k)[io;-e(k)] 



[iuj - e(k) - x A (k)S(k, iuj) 



(2) 



(3) 



Here %A(k) = 1 - l/{exp[(|e(k)| - A)/0.05A] + 1} cuts 
off contributions from e(k)| < A, and x'aO*-) = ^|A|X is 
nonzero only for |e(k)| near the cutoff A. Each diagram 
stands for two contributions since there are two ways to 
assign Ga and Sa to the internal lines. As in previous 
works [1-3] our calculation ignores higher order vertices 
and self-energy corrections (i.e Sa is set to zero). There- 
fore the RG flow has to be stopped before the interactions 
get too large and the lowest energy scales cannot be ac- 
cessed in a controlled way. 

We perform the RG numerically by dividing the Bril- 
louin zone into 32 patches. Each patch is centered on a 
"radial" line, as shown in Fig. (1). By standard Taylor 
expansion and power counting arguments, the most rel- 
evant vertex functions are the ones with incoming and 
outgoing ki on the fcrmi surface and the Wj = 0. In the 
following we only keep track of these most relevant vertex 
functions, which means we ignore the momentum depen- 
dence in directions transverse to the Fermi surface. For 
the instantaneous electron-electron scattering we thus 
approximate the value of VA(ki,k 2 ,k 3 ) by its value 
VA(k-F(i), kfO), ki?(7)), where kf(j') is the momentum 
on the fermi surface at the center of the jth patch. This 
leaves us with 32 3 couplings V(kp(i), kp(j), kp(Z)) to 
keep track of. We calculate the contributions to the 
flow of these couplings using the diagrams in Fig. (1). 
The frequency sums are performed analytically, and the 
loop momentum integration is performed by a sum over 
patches and integration along the radial lines. 





FIG. 1: a) One loop Feynman diagrams contributing to <9aVa- 
Spin is conserved along solid lines. The labels "I" and "R" 
represent the instantaneous and retarded interaction (the gray 
bars) respectively. While all interactions can be the "I" type, 
only those with both "I" and "R" labels can be retarded. A di- 
agram contributes to the flow of the instantaneous interaction 
when all interactions are the "I" type; otherwise it contributes 
to the flow of the retarded interaction, b) The 32-patch dis- 
cretization of the Fermi surface. Each patch (example in gray) 
is centered around one of the dashed lines. 



The retarded interactions have non-trivial frequency 
dependence. We assume them to have the form 



VR(ki,a;i;k2,W2;k3,a; 3 ) = - 



g(ki,k 3 )ff(k2,k 4 )£l_D 



(4) 

Here g(ki,k 3 ) is the electron-phonon matrix element 
for scattering an electron from momentum ki to k 3 , 
and Qjj is a characteristic phonon (Debye) frequency. 
The retarded interaction dies off for frequency transfer 
\u>i — w 3 | > flo- In this work we approximate this de- 
pendence by a step-function cutoff [8]. This approxima- 
tion throws away the full-frequency dependence of the 
retarded electron interaction used in other works[16] but 
makes the study of experimentally relevant fermi sur- 
face shapes calculationally tractable. This retarded in- 
teraction is represented by a new set of 32 3 "retarded" 
couplings, Vfi(ki?(i), k^O), k^(0)- Because the internal 
lines in each of the diagrams have frequency poles at or 
greater than the cutoff energy A while those of the ex- 
ternal legs are at zero energy, in the limit A >> Q,£> the 
leading contributions to the RG flow of the retarded in- 
teraction comes from diagrams in Fig. (1) where at least 
one interaction line is the "R" type. The contributions 
from diagrams with retarded interactions placed at ver- 
tices without an "R" label in Fig. (1) are down by a 
factor Qd/A. Because the phonon energies of interest 
in the cuprates are ~ 60 meV[17], in this work we in- 
tegrate the RG flow equations to a lower cutoff of 0.2t. 
Since the retarded interaction is present only when the 
frequency transfer is smaller than flo, diagrams with re- 
tarded interactions contribute only to the flow of retarded 
vertex function. Consequently the flow of instantaneous 
couplings is unaffected by the addition of retarded cou- 
plings, but the flow of retarded couplings is affected by 
the instantaneous interactions. 

The RG flow of the instantaneous quasiparticle scat- 
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FIG. 2: The renormalized instantaneous interaction 
VA=o.2t(kF(i), ki?(j), 1cf(0) plotted as a function of i and j for 
a fixed Z = 2. The solid vertical boxes (A) enclose couplings 
which promote commensurate and incommensurate magnetic 
ordering around (ir,ir). Examples of these scattering pro- 
cesses are shown on the right with each dashed arc enclos- 
ing a segment of the Fermi surface involved in the scattering 
processes. The dashed diagonal boxes (B) enclose couplings 
which promote d-wave superconductivity. The solid horizon- 
tal box (C) encloses couplings which lead to the Pomeranchuk 
instability. The circles (D) enclose couplings which lead to 4 
lattice constant charge-density-wave order. 



tcring alone with band parameter t" = has been ob- 
tained by Honcrkamp ct. al. We obtain similar results 
when including t" = — O.lt. Depending upon the doping 
level, several groups have found that the renormalization 
flow enhances instantaneous couplings which favor either 
SDW formation or d-wave superconductivity as the most 
prominent ordering tendency [1-3]. These can be seen in 
Fig. (2) where V\=o.2t 0^F{i), ^fO), is plotted as a 

function of i and j for a fixed 1 = 2. The couplings which 
favor SDW ordering have a positive amplitude and con- 
stant momentum transfer kp(Z)[l]. They show up 
as the vertical bands in boxes (A) of Fig. (2). The values 
k^O') — of these scattering processes are around 
(ir, it). In particular they include (n ± 27r/8a,7r), indi- 
cating "incommensurate" SDW ordering. The couplings 
which favor d-wave superconductivity are the diagonal 
lines enclosed in boxes (B) of Fig. (2). These scatterings 
occur in the cooper channel where kp(i) + = 0, 
and have an alternating sign structure indicative of d- 
wave symmetry. The fact that the couplings favoring 
superconducting pairing overlap with the couplings fa- 
voring SDW ordering signifies that physically these or- 
dering tendencies are linked. Since the phonon-mediated 
retarded couplings do not influence the flow of the instan- 
taneous interactions, the d-wave superconducting pairing 
seen above can not be due to the electron-phonon inter- 
action. In the present calculation, due to the choice of 
chemical potential and U, the spin-density waves are al- 
ways the leading ordering tendency. 

In addition to these two types of interactions, there 
also exist weaker scattering processes which favor other 
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FIG. 3: The three enhanced sets of retarded interactions, with 
their major types of scattering processes plotted as in Fig. (2). 
kF(Z) is fixed at Z = 2. Each set of retarded interactions has 
both positive and negative couplings. The sign structure is 
indicative of a x 2 — y 2 symmetry in momentum space. 



types of ordering. Couplings which lead to these order- 
ing tendencies are also indicated in Fig. (2). The first is 
the Pomeranchuk instability, which leads to a x 2 — y 2 - 
symmetry deformation of the Fermi surface [3]. The sec- 
ond leads to a charge-density wave order with ordering 
wavevector equal to the vectors connecting the parallel 
segments of the Fermi surface near the Brillouin zone 
face. These wavevectors correspond to a spatial period 
near four lattice constants. It is worth noting that the 
t' and t" terms in the band structure are responsible for 
the almost nested Fermi surface near the Brillouin zone 
face. They enhance the charge ordering tendency. 
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Now we investigate the flow of the retarded cou- 
plings. By experimenting with the symmetries of the 
initial scattering amplitude, we find three independent 
sets of retarded interactions that are most strongly en- 
hanced by the RG (Fig. (3)). In all three sets the scat- 
tering amplitude shows x 2 — y 2 symmetry in momen- 
tum space. This type of sign structure can arise from 
phonon mediation if g(k, k + Q) in Eq. (4) transforms 
as cos(k x ) — cos(ky). For the interaction in Figs. (3a), 
(3b), and (3c), Q - (±7r,±7r), Q ~ (±27r/4a, ±27r/4a), 
and Q ~ (±27r/4a, 0), (0, ±27r/4a), respectively. For ex- 
ample, it can arise from coupling to the half breathing or 
the Bi g [11] phonons. Note that the electron-phonon in- 
teraction with the Big buckling mode has no amplitude 
for momentum transfer (ir, 7r)[18] and so will not con- 
tribute appreciably to retarded interaction (3a). Consis- 
tent with previous results [12-1 5], the same analysis indi- 
cates that an s-symmetry electron-phonon coupling will 
be suppressed under RG. (Monte Carlo studies[15] find 
that s-symmetry electron-phonon coupling grows in the 
Hubbard model for large U > 6i; it would be of interest to 
see the behavior of x 2 — y 2 couplings.) Another common 
feature among the three sets is that the involved initial 
and final momenta are all within the antinodal region. 

We find that the most important contribution to the 
RG flow of the retarded interaction comes from the "bub- 
ble" diagram in the lower left of Fig. (l)(a), when one of 
vertices is an instantaneous interaction from the boxes 
(A) of Fig. (2) and the other is a retarded interaction. 
The x 2 — y 2 sign structure discussed above is essential 
for the growth of the scattering amplitude: the posi- 
tive(negative) retarded coupling combines with the in- 
stantaneous interactions to drive the negative (positive) 
retarded interaction to strong coupling. Since the boxes 
(A) enclose interactions that favor SDW ordering near 
(ir, it), we conclude that commensurate and incommen- 
surate magnetic fluctuations are responsible for driving 
the relevant retarded interactions to strong coupling. 

The key features of each set of retarded interactions in 
Fig. (3) appear as horizontal bands. In a horizontal band 
the momentum transfer kjr(i) — kjr(Z) is fixed, suggesting 
that these couplings promote some type of density wave 
order. However, the x 2 — y 2 sign structure makes the 
corresponding density wave order an unconventional one. 
The density wave order parameter driven by the retarded 
interactions in Fig. (3) has the form 

= ^/(k)<4 +Qa c k<7 >, (5) 

k,o- 

where /(k) transforms like cos(fc x ) — cos(k y ) under rota- 
tion. These are examples of generalized d-wave density 
orders, which in general can have bond averages < c\cj > 
with any complex phase, leading to both bond current 
and charge density orders [19]. We have verified via mean 
field calculations (details in a later publication) that all 
three interactions can promote real space modulations in 



both charge density and current. In particular, the in- 
teraction of Fig. (3c) can lead to periodic charge density 
modulations with period ~ 4a. Note that the chemi- 
cal potential /x can change the nesting wavevector and 
hence the period of modulations. Thus we find that in 
the cuprates the most important effect of the electron- 
phonon interaction is to promote density wave order. In- 
terestingly, for any order described by Eq. (5), the nodal 
quasiparticles will be unaffected. In particular the den- 
sity wave order in Fig. (3c) is consistent with the STM 
experiment of Hanaguri et al and the ARPES experi- 
ment by Shen et al. In addition, the fact that all the 
enhanced retarded interactions in Fig. (3) involve scat- 
tering processes with initial and final momentum states 
in the antinodal region is also consistent with the isotope- 
dependent ARPES study of Ref.[7]. 

Acknowledgement: We are in debt to J.C. Davis, K.P. 
McElroy, G.-H. Gweon, A. Lanzara, K. Shen, Z.-X. Shen, 
D. Scalapino, W. Hankc, and W. Metzner for many dis- 
cussions. DHL is supported by DOE grant DE-AC03- 
76SF00098. HCF was partially supported by a predoc- 
toral fellowship from the Advanced Light Source. 



[1] C. Honerkamp, M. Salmhofer, N. Furukawa, and T. M. 

Rice, Phys. Rev. B 63, 035109 (2001). 
[2] D. Zanchi and H. J. Schulz, Europhys. Lett. 44, 

235(1997); Phys. Rev. B 61, 13609 (2000). 
[3] C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 

(2000). 

[4] X. J. Zhou et al., Phys. Rev. Lett 92, 187001 (2004). 

[5] T. Hanaguri et al., Nature 430, 1001 (2004). 

[6] K. M. Shen et al., Science 307, 901 (2005). 

[7] G.-H. Gweon et al., Nature 430, 187 (2004). 

[8] G. T. Zimanyi, S. A. Kivelson, and A. Luther, Phys. Rev. 

Lett. 60, 2089 (1988). 
[9] A. Seidel, H.-H. Lin, and D.-H. Lee, Phys. Rev. B 71, 

220501 (R) (2005). 
[10] H. C. Fu, A. Seidel, H.-H. Lin, and D.-H. Lee (2006), to 

be published. 

[11] T. P. Devereaux, T. Cuk, Z.-X. Shen, and N. Nagaosa, 

Phys. Rev. Lett. 93, 117004 (2004). 
[12] J. H. Kim, K. Levin, R. Wentzcovitch, and A. Auerbach, 

Phys. Rev. B 44, 5148 (1991). 
[13] M. Grilli and C. Castellani, Phys. Rev. B 50, 16880 

(1994). 

[14] R. Zeyher and M. L. Kulic, Phys. Rev. B 53, 2850 (1996). 

[15] Z. B. Huang, W. Hanke, E. Arrigoni, and D. J. Scalapino, 
Phys. Rev. B 68, 220507(R) (2003). 

[16] S. W. Tsai, A. H. Castro Neto, R. Shankar, and D. K. 
Campbell, Phys. Rev. B 72, 054531 (2005). 

[17] T. Cuk, D. H. Lu, X. J. Zhou, Z. X. Shen, T. P. Dev- 
ereaux, and N. Nagaosa, Physica Stat. Sol. B-Basic Re- 
search 242, 11 (2005). 

[18] O. Jepsen, O. K. Andersen, I. Dasgupta, and S. Savrasov, 
J. Phys. Chem. Solids 59, 1718 (1998). 

[19] S. Chakravarty, R. B. Laughlin, D. K. Morr, and 
C. Nayak, Phys. Rev. B 63, 094503 (2001). 



